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Abstract 

We report on ab-initio calculations of the two-dimensional systems M0S2 and NbSe2, which 
recently were synthesized. We find that two-dimensional M0S2 is a semiconductor with a gap 
which is rather close to that of the three dimensional analogue, and that NbSe2 is a metal, which 
is similar to the three dimensional analogue of this compound. We further computed the electronic 
structure of the two-dimensional hexagonal (graphene like) lattices of Si and Ge, and compare them 
with the electronic structure of graphene. It is found that the properties related to the Dirac cone 
do not appear in the case of two-dimensional hexagonal germanium, which is metallic, contrary 
to two-dimensional hexagonal silicon, which has an electronic structure very similar to the one of 
graphene, making them possibly equivalent. 

PACS numbers: 71.15.Mb, 71.20.-b 



1 



I. INTRODUCTION 



During several years, it was believed that the existence of free-standing two dimensional 
crystals was impossible, because they would be unstable and ultimately return to a three 
dimensional object^. However, it appeared recently that isolated sheets of graphene could 
be obtained by mechanical exfoliation of a graphite crystal, which therefore proves the pre- 
dictions of Rcf 1 2 to be inaccurate 3 . The experimental accomplishment of synthesizing 
two-dimensional crystals has lead to the emergence of a truly new physics, since the particu- 
lar properties of graphene can be considered as a bridge between quantum electrodynamics 
(QED) and condensed matter physics^. Indeed, near the K point of the Brillouin zone 
(BZ), the one-particle energy dispersion is linear with the momentum, and therefore the 
corresponding quasiparticles could be described by a Dirac-like Hamiltonian. Then, QED 
properties can be studied by investigating the electronic structure of graphene. 

Later on, it was shown^ that the same technique could be used to obtain other com- 
pounds, opening the path to the investigation of a large number of two dimensional crystals. 
However, it seems that most of the experimental and theoretical efforts in this area are 
still focusing exclusively on graphene, primarily because of the high quality of samples that 
could be obtained. In the present paper, we investigate by means of ab-initio calculations 
the electronic structure of M0S2 and NbSe2, as they have been synthesized experimentally^, 
as well as hypothetical two dimensional hexagonal crystals (with a graphene like structure) 
made of Si and Ge. 



II. COMPUTATIONAL DETAILS 

To perform the calculations, we have used density functional theory (DFT)&£ as imple- 
mented in the code VASP (Vienna Ab-initio simulation package }2* 9 , within the framework 
of the PAW (projector augmented waves) method 10 . The Local Density Approximation^ 
(LDA) as well as the Perdew Burke Ernzerhol™ variant of the generalized gradient approx- 
imation (GGA) were used for the exchange-correlation potential. A cut-off of 500 eV was 
used for the plane-wave expansion of the wave function to converge the relevant quantities. 
For Brillouin zone integrations, a mesh of 40 x 40 x 1 fc-pointsi 3 - was used for M0S2 and 
NbSe2, while a mesh of 20 x 20 x 3 was sufficient to describe two-dimensional silicon and 
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germanium. 

Concerning the crystal structure M0S2 and NbSe2 we used the experimental bulk values 
since it was noticed in Ref. ||J that 2D crystals in these cases remains very close to their 3D 
parents. In this case, GGA was used. For the two dimensional hexagonal crystals made of 
Si and Ge, we optimized the lattice parameter using either the LDA or the GGA. Moreover, 
it appeared that the band structure obtained using either LDA or GGA were very similar, 
so we have chosen to show only results corresponding to LDA. For all studies here, the 
'c' parameter was taken to be large enough to ensure that no interaction remains between 
layers, making them effectively isolated 2D objects. 



III. THE ELECTRONIC STRUCTURE OF TWO-DIMENSIONAL MOS 2 AND 
NBSE 2 



A. MoS 2 

There is a big interest in transition metal dichalcogenides, since due to their layered 
structure, they have extremely anisotropic properties, and therefore an intercalation process 
is easy to conduct. In particular, M0S2 is used in the technology of Li batteries. To 
investigate purely two-dimensional M0S2, we have used the bulk lattice parameter (a = 3.16 



A), as reported in Ref Il4j. Since reduced dimensionality can sometimes leads to magnetic 
behavior in systems which are not magnetic in bulk, we checked the possibility to have 
a spin-polarized ground-state, but it was found that two-dimensional M0S2 remains in a 
non-magnetic state. 

Our calculated DOS is presented in Fig [U As for three-dimensional M0S2, the two- 
dimensional variant of this compound is semiconducting: the bands on each side of the 
band gap are derived mainly from the Mo-d states, which is what a calculation of three 
dimensional M0S2 also shows 15 . Also, a large gap exists in the occupied states, separating 
the S-s states (between -14 and -12 eV) and the hybridized Mo and S-p states (between 
-6 and eV). In Fig. [21 we present the energy bands of two-dimensional M0S2. We see 
again a lot of similarities with three dimensional M0S2: the bands around the band-gap are 
relatively flat, as expected from the 'd' character of the electron states at these energies. 
The band gap, which has a value of 1.78 eV, is direct, and occurs at the high-symmetry 
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FIG. 1: The density of states of two-dimensional M0S2. The Fermi level is put at zero eV. 




FIG. 2: The band structure of two-dimensional MoS2- The Fermi level is put at zero eV. 
point K, whereas for three dimensional M0S2 the band gap is indirect^. 



B. NbSe 2 



As for MoS2, the NbSe2 structure is strongly layered, with each Nb layer being sandwiched 
between two Se layers, where weak van der Waals forces are holding the whole geometry 
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together. NbSe2 is known to be a prototype to study charge density waves (CDW), and is 
also a superconductor. Here we focus only on the electronic structure of two-dimensional 
NbSe 2 and a lattice constant of 3.45 A— was used to perform the calculation. 
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FIG. 3: The density of states of two-dimensional NbSe2- The Fermi level is put at zero eV. 

As for M0S2, we have checked the existence of a magnetic solution, since the non spin- 
polarized DOS (Fig. [3]) is very peaked at the Fermi level, but we found again a non-magnetic 
ground state. However, contrary to M0S2, two dimensional NbSe2 is metallic, just like the 
three dimensional version of this compound. The DOS of two-dimensional NbSe2 shows a 
lot of similarities with the one of M0S2, with a large hybridization between Nb and Se states. 
In addition, the band state which is pinned at the Fermi level is derived primarily from Nd 
d-orbitals, and this band state becomes separated in energy from all other states (Fig. H]). 

IV. THE ELECTRONIC STRUCTURE OF TWO-DIMENSIONAL SI AND GE 

With four 'sp' electrons in the valence band, silicon and germanium (possibly together 
with BN) are probably the closest to carbon from a chemical point of view. Therefore, it is 
of interest to see if in their two-dimensional form they could present some similarities with 
graphene, and we have investigated this possibility here. 

In Table HJ we present our computed lattice parameters (see computational section) for 
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FIG. 4: The band structure of two-dimensional NbSe2- The Fermi level is put at zero eV. 



System 


LDA 


GGA 


Si 


3.860 A 


3.901 A 


Ge 


4.034 A 


4.126 A 



TABLE I: Lattice parameters of hexagonal Si and Ge computed with either LDA or GGA. 

hexagonal Si and Ge. As expected, they are larger than the one of graphene (2.46 A) because 
of the larger radius of Si and Ge in comparison with C. Also, because of the well-known 
overbinding of LDA, GGA gives lattice parameters which are larger than the ones given by 
LDA. 

We used these lattice parameters for the calculation of the corresponding densities of 
states and band structures. As seen in Figs [5] and El two-dimensional silicon shows a lot 
of similarities with graphene: in particular the gap is also closing at the K point of the 
Brillouin zone, and the dispersion around this point is linear. In contrast, two-dimensional 
germanium is quantitatively different from graphene: a conduction band is partially filled by 
electrons around the F point and therefore the bands at the K point are shifted up in energy, 
so that two-dimensional germanium is not a zero-gap material, but rather a poor metal. For 
three dimensional germanium, it is known that LDA/GGA gives a non band-gap behavior, 
and a more advanced technique like the GW approximation^&^ is required to obtain a 
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FIG. 5: The density of states of two-dimensional Si. The Fermi level is put at zero eV. 




FIG. 6: The band structure of two-dimensional Si. The Fermi level is put at zero eV. 



band gap. However, our calculations with DFT show that two-dimensional germanium is 
much more metallic than three-dimensional germanium, and therefore the metallic nature 
is here much more robust and is likely to be preserved even if a quasiparticle theory (like 
the GW approximation) would be used. 
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FIG. 7: The density of states of two-dimensional Ge. The Fermi level is put at zero eV. 




FIG. 8: The band structure of two-dimensional Ge. The Fermi level is put at zero eV. 
V. CONCLUSION 

We have studied by means of ab-initio calculations the electronic structure of two- 
dimensional MoS 2 and NbSe 2 , as well as hypothetical graphene like structures of Si and 
Ge. We have found that two-dimensional silicon might from an electronic structure point of 
view be equivalent to graphene, with a linear dispersion of the electronic structure around 
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the K-point. The possible advantage with Si in this regard is that it probably is more easily 
interfaced with existing electronic devices and technologies. The obvious disadvantage is 
that sp 2 bonded Si is much less common than for C, and the synthesis of Si in a graphene 
like structure is extremely demanding and is likely to represent a meta-stable material. The 
electronic structure of Ge in the graphene structure results in a metallic behavior, and the 
electronic structure of two-dimensional M0S2 and NbSe2 is quite similar to that of the three 
dimensional counterparts, with M0S2 being a gapped system with a direct gap and NbSe2 
a metal. 
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